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1.  Introduction 

y  juC. 

JJ*  shall  be  interested  here  .in  finite  element  discretizations  of 

i  /\ 

problems  involving  an  incompressibility  condition.  As  model  problems  we 
consider  the  Stokes  equations  for  the  flow  of  a  viscous,  incompressible  fluid 
and  the  equations  of  linear  plane-strain  elasticity  for  the  deformation  of 
an  isotropic,  nearly  incompressible  solid.  In  both  cases  the  incompressibility 
condition  takes  the  form  of  a  divergence  constraint.  Although  this  is  the 
most  simple  formulation,  the  proper  understanding  of  how  an  approximate 
method  satisfies  the  constraint  represents  an  important  step  towards  the 
understanding  of  more  complicated  situations,  involving  e.g.  the  Navier- 
Stokes  equations  or  the  equations  of  nonlinear  elasticity.  The  finite 
element  methods  we  study  have  the  property  that  the  approximations  to  the 
velocities,  respectively  to  the  displacements,  are  continuous;  such  methods 
are  generally  referred  to  as  conforming. 
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2.  The  Stokes  equations 


Let  ft  be  a  bounded  polygonal  domain  in  the  plane,  and  let  U  and 


P  solve 


(2.1) 


-AU  +  VP  *  F  in  ft  , 

7*U  »  0  in  ft  ,  and 


U  *  0  on  3ft  . 


Here  £  =*  (U^.U^)  represents  the  fluid  velocities  and  P  the  pressure;  the 
viscosity  has  been  set  to  1  .  To  simplify  the  exposition  we  are  assuming 
homogeneous  boundary  data  on  3ft  .  For  the  linear  problem  (2.1)  we  can 
convert  inhomogeneous  boundary  data  into  an  external  force  term  £  .  However, 
for  a  nonlinear  problem  one  must  deal  directly  with  the  inhomogeneous 
boundary  data,  cf.  Gunzburger  &  Peterson  [17].  Other  -  ble  boundary 
conditions  could  involve  the  normal  fluid  stresses 

2  f3U  3U 

-P"i  >  1*1-2’ 

but  since  stress  boundary  conditions  are  physically  much  more  frequent  when 
dealing  with  solids,  we  shall  reserve  these  for  our  formulation  of  the  boundary 
value  problem  for  the  equations  of  elasticity.  For  regularity  results  con¬ 
cerning  the  solution  to  (2.1)  on  a  polygonal  domain  see  Kellogg  &  Osborn  [20]  and 
Osborn  [29],  Note  that  P  is  determined  only  up  to  an  additive  constant. 


w **  wj TAru.jmri*m.T**?.w't? •? *.vvav  ^  tv-v-v^.v: 
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3.  The  equations  of  linear  elasticity 

As  before  let  Q  be  a  bounded  polygonal  domain  In  the  plane  and 
consider  the  problem 


(3.1) 


I 


-AU  - 
—  l-2v 

2 

l  - 


7(V«U) 


F  in  A 


j-1 


iJ®>nJ  +  I^77^ni 


and 

on  90  , 


1-1.2  , 


ii  *  (n^,^)  here  denotes  the  outward  unit  normal  to  0  ,  and  e^(lJ)  is 
the  usual  symmetric  strain  tensor 


Elj(£) 


1 

2 


The  equations  (3.1)  are  the  equations  of  isotropic,  plane-strain  linear  elasticity 
corresponding  to  the  domain  0.  0<v<l/2  is  a  material-dependent  con¬ 
stant,  the  so-called  Poisson's  ratio,  which  describes  the  compressibility. 

The  other  constant,  the  shear  modulus,  that  is  needed  in  order  to  characterize  fully 
an  isotropic  material  has  been  absorbed  into  the  external  load  F  and  boundary 
load  £  .  If  U  is  a  solution  to  (3.1)  then  the  vector 


U  -  (U^.O) 


solves  the  equations  of  3-dimensional  isotropic  elasticity  on  the  domain 

n  x  e  , 

with  external  load  £  *  (F,0)  ,  independent  of  ,  and  boundary  load 
fL  *  (fiL»0)  »  independent  of  x^  ,  on  the  vertical  boundary  9ft  x  ]R  .  This 
is  the  reason  for  the  notion  of  "plane-strain"  In  connection  with  (3.1). 
Values  of  v  near  1/2  correspond  to  a  nearly  incompressible  material. 
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Both  the  problem  (2.1)  and  the  problem  (3.1)  can  (on  a  simply  connected 

domain)  be  reduced  to  the  solution  of  the  biharmonic  equation  with 

Dirichlet  boundary  conditions.  For  (2.1)  this  follows  through  the 

introduction  of  a  stream  function  such  that  U  =  7  x  y  ;  the  corresponding 

boundary  conditions  for  y  are  homogeneous  Dirichlet.  For  (3.1)  one  first 

subtracts  a  particular  solution  corresponding  to  the  external  load  £  so 

that  the  resulting  system  has  a  vanishing  external  load.  For  such  a  system 

2 

one  may  introduce  an  Airy  stress  function  $  (  A  $  =  0)  so  that 


2  2 

r  ^  > 

3 

a  ^  a  —  _ 

3 

3x^ 

4 

CT22  »  3X;l3x2$  a12  ’ 

3x_ 

where 


ij 


E 

1+v 


(U)+<5 


7*U 


ij  -  ij  l-2v 


denote  the  stresses  corresponding  to  the  displacement  IJ  (E  is  the  so-called  Youngs 
modulus)  .  The  boundary  conditions  for  $  are  in  this  case  (inhomogeneous)  Dirichlet. 

Regularity  results  and  a  priori  estimates  for  solutions  to  (3.1)  (as 
well  as  (2.1))  may  thus  be  derived  from  the  properties  of  solutions  to  the 
biharmonic  equation  on  a  polygonal  domain  (cf.  Grisvard  [15]). 


4.  Variational  formulations 

In  order  to  introduce  finite  element  discretizations  of  the  equations 

(2.1)  and  (3.1)  we  have  to  cast  them  in  a  variational  (or  weak)  form.  This 
is  done  following  the  standard  3-step  recipe: 

1)  multiply  each  differential  equation  by  a  suitable 
test  function 

2)  integrate  the  result  over  ft 

3)  integrate  by  parts  (to  taste). 

By  multiplication  of  the  first  equation  in  (2.1)  by  v  (vanishing  on  3ft)  the 
above  three  steps  lead  to 

a(U,v)  +  b(v,P)  =  (F,v) 

2 

where  (F,v)  is  the  usual  [I^ft)]  inner  product, 

(4.1)  a(U,v)  *  2  l  ,  and 

•'ft  i,j  1J  J 

* 

(4.2)  b(v,P)  =  -  V*vPdx  . 

'ft 

Multiplying  the  second  equation  in  (2.1)  by  a  suitable  function  q  and 
integrating  over  ft  ,  we  get 

b(U,q)  =  0  . 

The  appropriate  spaces  of  "test"  functions  v  and  q  are  given  by 

(4.3)  H^ft)  x  h1^)  and  L2(ft)  . 

°1 

Here,  H  (ft)  is  the  standard  Sobolev  space  of  functions  whose  gradients  are 
square  integrable  and  whose  traces  vanidi  on  3ft  .  With  these  spaces  our 
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"variational  form"  of  (2.1)  is 


(4.4) 


Find  U  €  [H1^)]2  and  P  €  Lo(0)  such  that 


a(U,v)  +  b(v,P)  =  (F, v)  for  all  v  €  [h\o)]2 


b(U,q)  *  0  for  all  q  6  . 


In  order  to  find  the  variational  formulation  of  (3.1)  we  multiply  the  first 

1  2 

equation  by  v  €  [H  (0) 1  (not  vanishing  on  80  )  and  integrate  by  parts. 
Because  of  the  form  of  the  boundary  conditions  this  leads  to 


a(U,v)  +  b(v,  -  Y32\7  7‘-^  =  *— ■ ^  +  2<E.*0?  » 


where  <•,•>  is  the  [L^OO)]  inner  product.  The  differential  equation  (3.1) 
therefore  has  the  weak  form 


(4.5) 


Find  U  €  [H^(0) ]2  such  that 


a(U,v)  +  b(v,  -  j- 2^  V*U)  =  (L>Z)  +  2<£»Z> 


This  may  be  rewritten  as 


Find  IJ  €  [H^(0)]2  and  P  €  1^(0)  such  that 


a(£,v)  +  b(v,P)  *  (£,v)  +  2<&,v>  for  all  v  €  [H^"(0)]^ 


^ -  (P,q)  +  b(U,q)  =  0 


for  all  q  €  I^ft) 


Setting  v  ■  1/2  this  gives  "formally"  a  Stokes  equation  of  the  form  (4.4) 
only,  of  course,  with  different  boundary  conditons.  Indeed  the  equations  of 


elasticity  for  v  near  1/2  may  be  viewed  as  a  penalized  version  of  a 
Stokes  problem  (for  more  details,  see  e.g.  Temam  [35]).  Note  that  the 
equation  (4.5)  only  has  a  solution  provided  the  loads  £  and  £  are 
statically  admissible,  i.e.  provided 

(F,R)  +  2<g,R>  -  0 

for  any  rigid  motion  R  .  The  solution  _U  is  also  only  determined  modulo  a  rigid 
motion;  whenever  we  discuss  the  problem  (4.5)  we  shall  avoid  this  non-uniqueness 
by  thinking  of  functions  as  equivalence  classes  modulo  rigid  motions 
(remember,  a  rigid  motion  is  one  for  which  (R)  =  0  ,  or  ^.(x^.x^)  =  (-yx2+a, 
YXj+8)  .  In  this  paper  we  shall  always  assume  that  £2  is  connected  and  for 
the  Stokes  equations  we  make  P  unique  by  imposing 

» 

P  d  x  =  0  . 

o 

This  replaces  L„(ft)  in  (4.4)  by  L„(Q)  =  L9(fi)  (1  {q  :  q  dx  =  0) 


5.  Stability  of  finite  element  approximations 

Given  a  variational  formulations  such  as  (4.4)  or  (4.5)  the  finite 
element  method  consists  of  choosing  finite  dimensional  subspaces 


^SIH^Q)]2  ,  nhcL2(Q) 


V.  c  [HX(fi)]2  , 


and  replacing  everywhere  the  infinite  dimensional  spaces  by  their  finite 

*  0 

dimensional  counterparts.  The  spaces  and  are  typical  ^ade  up  of 

piecewise  polynomial  functions  on  some  triangulation  ^  ;  h  denoti  :he 
mesh  size.  For  the  Stokes  problem  the  discrete  version  thus  becomes 


Find  U,  €  V  and  P  €  II,  such  that 
— n  h  h  h 


(5.1) 


a (U^ , v )  +  b(v,Ph)  «(F,v)  for  all  v  €  Vh 


bO^.q)  -0  for  all  q  €  JI 


and  for  the  equations  of  elasticity  (4.5)  it  reads 


Find  i  such  that 


(5.2) 


a (U^ » v)  +  b(v, -  -_2~  7.^)  -  (F,v)  +  2<£,v> 


for  all  v  €  V,  . 
—  n 


*  12 
Indeed  shoi  be  a  subset  of  [H  (ft)]  /{Rigid  motion^  out  this  is  not 

explicitly  menti  :  our  convention  is  to  identify  a  func  .on  with  ita 

equivalence  clas  iulo  rigid  motions,  whenever  appropriate. 


v  *.'.v  s'  >\v  ■ 


As  an  example  we  consider  a  uniform  triangulation  which  locally  is  as  shown  in 
Fig.  1.  Perhaps  the  simplest  finite  dimensional  subspaces  are 


0  o 

[C  piecewise  linear  functions]^  ,  or 


°  0  2 
(5.3)  =  [C  piecewise  linear  functions  that  vanish  on  3ft] 

and  correspondingly 

II  =  piecewise  constants  having  integral  zero  (over  ft). 


One  reason  why  finite  element  methods  for  the  equations  of  elasticity 
(near  incompressibility),  or  the  Stokes  equations,  are  intriguing  is  that  the 
choices  of  spaces  (5.3)  d£  not  work  for  the  mesh  in  Fig .  1 !  In  the  case  of 
the  discrete  equations  of  elasticity  (5.2)  the  relative  error  as  v  -*■  1/2  for 
fixed  h  approaches  a  constant,  which  is  bounded  away  from  0  independently 
of  h  .  For  the  discrete  Stokes  equation  (5.1),  =  0  for  any  h  .  Below 

we  give  an  explanation  of  this  phenomenon  in  the  case  of  the  Stokes  equations. 
Since  the  Stokes  equations  are  the  "limit"  as  v  -*•  1/2  of  the  equations  of 
elasticity  this  also  intuitively  explains  the  lack  of  uniformity  in  v  of 
the  accuracy  of  the  approximation  (5.2)  to  the  equations  of  elasticity.  (Note, 
however,  the  difference  in  our  boundary  conditions.) 

The  second  equation  in  (5.1)  requires  that  lies  in  the  subspace 


(5.4) 


Zh  =  ^  v v,  :  b(v_,  q)  =  0  Vq  €  II,  } 


and  part  of  (5.1)  may  thus  be  restated 


(5.5) 


Find  U,  €  Z,  such  that 
-h  h 


a(U^,v)  =  (J\v)  for  all  v  6  Z^ 


The  reason  for  the  deficiency  of  the  choice  of  spaces  and  IT,  ,  given  by 


i 


(5.3)  is  that  the  corresponding  Z^  on  the  mesh  £  consists  of  the  zero 
element  only.  This  can  be  seen  as  follows:  we  first  observe  that 


Zh  =  V  *  {C  piecewise  quadratics  all  of  whose  first 
derivatives  vanish  on  9ft}  , 


with  the  curl  operator  V  x  given  by 


A  C  piecewise  quadratic  whoose  first  derivatives  vanish  on  9 ft  can  be 

2 

extended  by  constants  onto  R  \ft  .  By  subtracting  the  constant  which  is  attained 

2  1 

in  the  unbounded  component  of  R  \  ft  we  thus  obtain  a  C  piecewise  quadratic 
with  compact  support.  On  the  mesh  shown  in  Fig.  1  there  is  only  one  piecewise 
quadratic  with  compact  support  (cf.  Morgan  &  Scott  [26],  Chui  &  Wang  [ 7  ]) 
namely  the  constant  0!  Consequently  it  follows  that 


Z  =  V  *  (constant  function} 
n 

=  {0}  . 


The  formulation  (5.5)  mimics  the  following  equivalent  version  of  the  Stokes 
problem 


(5.6) 


Find  U  ?  Z  such  that 


a(lJ,v)  =  (F,v)  all  v  €  Z 


where  Z  denotes  the  subspace 


Z  ■  (v  €  [H^(ft)]^  ;  b(v,q)  ■  0  Vq  €  I^fi)} 


(v  €  [H^(ft)]^  :  7»v  *  0} 


12 


Note  that  (5.6)  represents  the  well  known  Hodge  decomposition  (cf.  Temam  [35]). 
The  space  Z^  ,  as  defined  in  (5.4),  is  of  course  not  necessarily  contained  in 

O 

Z  ,  but  in  many  interesting  cases  this  inclusion  holds,  i  ■  e .  Z^  =*  fl  Z  . 
Since  a(*,*)  is  a  positive  definite  symmetric  form  that  coerces  the  H^-norm 

O  1 

(on  H  (fl))  it  follows  in  this  case  that 


|U-IL  I!  <  C  inf  llu-z.ll 

H  (fl)  z€Z.  H  (fl) 

—  n 


We  would  thus  obtain  a  quasioptimal  velocity  approximation,  ,  if  the  follow¬ 
ing  condition  were  to  hold: 


(5.7) 


For  any  U  €  Z  , 


inf  ||  U-z||  i  S  C  igf  ||U-v|| 

z€Z,  H  (fl)  v€V,  H  (fl) 

—  h  —  h 


with  C  independent  of  U  and  h  . 


Definition  I 


A  family  of  closed  (not  necessarily  finite  dimensional)  subspaces 
1  2 

Wh  c  [H  (fl)]  is  called  divergence-stable  if 


i)  the  spaces  V*W^  are  closed  in  L^tfl)  , 


ii)  3c  >  0  ,  independent  of  h  ,  such  that 


sup  —  L 

w€W.\{0}  ||w|| 

h  H^fl) 


*  C  II  q  I 


Mfl) 


for  all  q  €  V*W  . 

n 


2 


Remark  5.1 


The  definition  of  divergence-stability  as  given  above  is  equivalent  to 
the  requirement  that  there  exists  a  uniformly  bounded  maximal  right  inverse 
for  the  divergence  operator  on  the  spaces  ,  i.e.  there  exists  a  family 
of  linear  operators 

L  :  V*W,  -*■  W, 

.  _  n  h  h 

such  that 

i)  V*(Lhq)  -  q  ,  *q  €  V*Wh  and 

ii)  MLhqNi,n  *  c^q^o,a  ’ 

with  a  constant  C  that  is 
independent  of  h  and  q  . 

The  condition  (5.7)  and  the  concept  of  divergence-stability  are 
intimately  related  as  shown  by  the  following  result. 

Proposition  5.1 

The  spaces  fl  Z  and  satisfy  (5.7)  for  any  £  €  Z  ,  with 

O 

constant  C  that  is  independent  of  h  and  U.  ,  if_  and  only  if  is 

divergence-stable . 

Proof: 

o  o 

Assume  that  Z.  *  V,  fl  Z  and  V,  satisfy  (5.7).  For  any 
n  n  n 

o 

q  €  V*V^  there  exists 

v  €  [H1^)]2  , 


such  that 


V«v  -  q  and  ||vj|  .  5  C||q|| 

h  (n)  l2(«) 


(cf.  Temam  [35]). 


Since  q  €  V*V,  there  also  exists  v.  €  such  that 
h  — h  h 


V*v^  ■  q 


Hence 


v-  ^  6  Z  , 


and  therefore  there  exists  z,  €  Z,  with 

h 


(5.8) 


i  $  c  f2f  Hr-Zh-wll  i 

H  (0)  w€$h  H  (fi) 


*  C  llv.ll 

HA(£2) 


The  estimate  (5.8)  immediately  leads  to 


(5.9) 


1  *  C  Hill  l 

^  H  (Q)  Hx(n) 


*  c  ||  q  | 


l2(«)' 


Since  also 


(5.10) 


7-(W  *  1 


the  field 


-  (v,  +  z,)  €  V,  satisfies 
— h  — h  h 


(5.11) 


iV  l 
h  <n) 


HqllL2(n) 


^  c  ||  q  | 


iHhll  i 
^  h  (n) 


o 


(0)  • 


As  a  consequence  of  (5.9)-(5.11)  we  conclude  that  V  is  divergence-stable 


The  proof  in  the  opposite  direction  is  a  simple  consequence  of  the 
results  of  Babuska  [  2  ]  and  Brezzi  [  6  ];  for  completeness  we  include 

O  • 

the  details  of  a  proof  here.  For  each  v  €  ,  ^et  ^h  *  Vh  Z  be  itS 

orthogonal  projection  with  respect  to  the  inner-product  a(*,*)  ,  and  define 


qh  €  V.Vh  via 


b(w,q,  )“  a(v-z  ,w)  for  all  w  €  V 


9 

Since  V.  is  divergence-stable,  q^  is  well  defined  and  satisfies 


|UI  S  cfc-sjl 

^  L  (0)  ^  H  (SI) 


Now  by  definition  of  z_  an<^  ^  » 


c  WsrS^W  i  *  a<v-z^,v-z^) 

^  H  (ft) 

-  a(v-z,  ,v) 


b(v,qft) 


»  b(v-U,qh) 

*  C  ||  v-Uj|  ||qhi| 


H  (ft)  L2(fi) 


S  C  ||  v-U.ll  1  lt-*JI  l 

hV)  Ha(B) 


for  any  }J_  £  Z  .  Thus  we  have,  for  any  V  €  Z  , 

||v-zh[|  !  ?  C||v-U||  1  ^ 

IT(ft)  H  (ft) 


By  the  triangle  inequality,  this  means  that,  for  U  €  Z  and  v  €  , 


provided  is  the  projection  (relative  to  the  inner  product  a)  of 

O 

v  onto  0  Z  .  The  estimate  (5.12)  clearly  implies  (5.7).  (Note  that  at 

O 

no  point  in  this  proof  did  we  use  that  is  finite  dimensional). 


Remark  5.2 

Proposition  5.1  does  not  assert  that  divergence-stability  is  alwaya 
necessary  in  order  to  get  optimal-rate  velocity  approximations.  Indeed  optimal- 
rate  velocity  approximations  are  in  certain  cases  achieved  with  piecewise  polynomial 

O 

spaces  that  are  not  divergence-stable  (ef.  Remark  6.1). 


We  now  suppose  that  €  Z^  is  known  through  the  solution  of  (5.5),  and 
we  consider  determining  €  JI^  from  the  first  equation  in  (5.1) 

O 

Cj.IJ)  l»(v,l*h)  =  (l-'.v)  -  .i(llh,v)  for  .ill  y  '  V  . 


In  order  for  P^  to  be  unique  it  is  necessary  and  sufficient  that  there  exists 
c  >  0  such  that 


(5.  14) 


l<  ( v ,  >1 ) 

sup  - : - 

!MI  , 

H  (Si) 


I.,(Q> 


for  all  q  €  .  If  we  furthermore  want  ||PjJ|L  ^  uniformly  bounded  by 

C  ( 11 P || .  +  Hull  i  )  independently  of  h  then  if  suffices  to 

L2U  ;  H  (Si) 

require  that  the  constant  c  in  (5.14)  is  independent  of  h  . 

We  note  that  although  it  is  mathematically  simple  to  describe  it  is  in 


practice  not  always  obvious  how  to  use  the  equations  (5.5)  ami  (5.13)  to  solve 
for  i.  and  P.  (cf.  Glowinski  [14]). 


a 


r-\'. 


& 


$ 


r$ 


a 


R 
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Definition  II 


1  2 

Let  W.  and  Q,  be  two  families  of  closed  subspaces  of  [H  (Q)] 
n  a 


and  respectively.  Tne  family  is  said  to  be  divergence-stable 


relative  to  if  3c  >  0  ,  independent  of  h  ,  such  that 


L>(w,q) 


sup 


w(W  \  {0}  1 1  w  1 1  . 

H  (fi) 


:  •:  II  q  I 


L0(Q) 


for  all  q  €  Qh  . 
Remark  5.3 


By  comparison  of  the  definitions  I  and  II  we  note  that  a  family  of 


closed  subspaces  c  [H  (H)]  is  divergence-stable  if  and  only  if  it  is 


divergence-stable  relative  to  the  family  7*  . 


Using  the  concepts  of  divergence-stability  introduced  in  definitions  I 
and  II  the  question  of  stability  of  finite  element  approximations  to  the 
Stokes  equations  (2.1)  and  the  equations  of  elasticity  (3.1)  is  well  understood. 
Proposition  5.2 


o  7 

Suppose  that  the  family  of  spaces  V  £  [H  (Q) ]“  is  divergence-stable 


L,(0)  .  Let  (U.P)  €  fH1]2  *  4  and  (U^)  e  \  x  \ 


relative  to  the  family  14  c 
Jointly  the  sol  tit.  ions  to  ('*./,)  and  ('’.!)  respect  IvoLv;  then 


llu-u.ll  L  +||p-pJI 

H  (ft) 


l2(o) 


$  C  (  inf  I! U-v.ll  ,  +  inf  || P-qlL 


v€V, 
—  n 


h  (a)  q€n. 


and  C  is  independent  of  h  (and  F) . 

Proposition  5.2  follows  directly  from  the  results  of  Babulka  [2]  or 
Brezzi  [6]«  It  reduces  the  question  of  convergence  to  one  of  approximation  theory 
and  this  inherently  leads  to  contradictory  requirements:  on  one 


hand  the  spaces  should  have  good  approximation  properties  -  on  the  other 

hand  they  should  not  be  too  big,  since  that  would  jeopardize  stability.  One 

O 

very  natural  choice  of  II,  seems  to  be  IT,  «  V*VL  . 

h  h  h 

Proposition  5. 3 

o  °1  2 

Suppose  that  the  family  of  spaces  £  [H  (SI)]  is  divergence- 
stable,  and  choose  n  =  V-V  .  Let  (U,P)  €  [ H1  ] 2  -  L0  and  (U.  ,P  )  €  v.  *  il 


*  h  v  vh  h 


be  the  solutions  to  (4.4)  and  (5.1)  respectively;  then 


i>  l|U-£hll  i  *  C  inf  ||U-v[|  and 

H  (SI)  v€Vh  H  (SI) 


ii)  llP-P  ||  <  C  inf  ||U-v||  +  inf  ||P-qJ| 

L2(SI)  \^v6Vh  n  (SI)  q€nh  4  W 


and  C  is  independent  of  h  (and  F) . 

Proof 

This  follows  immediately  from  Prop.  5.1,  Prop.  5.2  and  the  analysis 
leading  to  the  definition  I. 

In  the  following  we  shall  denote  by  €  [H  L]2  *  ^  Vh  the  solutions 
to  (4.5)  and  (5.2),  thus  emphasizing  the  dependence  on  0  <  v  <  1/2  of  the 
solution  to  the  problem  of  elasticity. 

Proposition  5.4 

1  2 

Suppose  that  the  family  of  spaces  £  [H  (SI)  ]  is  divergence- 
stable.  Then 


H  (SI) 


*  C(  inf  ||UV-v||  +  inf  H73-  ?’UV-q(|  . 

\v€vh  n  (si)  q€S7.Vh  2U; 


with  a  constant  C  that  is  independent  of  h  and  v  (and  F  and  g) . 


Proof 
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For  a  fixed  0  <  v  <  1/2  this  estimate  with  a  constant  C 
possibly  depending  on  v  ,  follows  directly  from  the  results  of  Babu§ka  [2  ]. 
The  fact  that  C  may  be  chosen  independently  of  v  ,  provided  the  family  V 

h 

is  divergence-stable,  is  e.g.  proven  in  Vogelius  [36], 


Several  ways  to  remedy  the  deficiency  of  the  spaces  in  (5.3)  have  been 
suggested  and  analyzed  in  the  literature.  A  non-conforming,  piecewise  linear, 
triangular  element  was  proposed  by  Crouzeix  and  Raviart  [10]  which,  together 
with  the  1^  in  (5.3)  yields  optimal  approximation  (second  order).  For 
computational  experience  with  this  element,  see  Pritchard,  Renardy  &  Scott  [31] 
and  references  therein.  Increasing  can  also  be  achieved  by  choosing  a 

special  triangulation,  as  shown  in  Figure  2,  although  this  Is  by  no  means 
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It  was  discovered  by  Powell  [30]  that  the  space  of  piecewise 

quadratics  on  the  mesh  in  Figure  2  has  good  approximation  properties  -  there 
is  roughly  one  degree  of  freedom  per  singular  vertex.  One  must  reduce 
in  (5.3)  appropriately  at  singular  vertices,  cf.  later  in  this  paper.  For 
computational  experience  with,  and  analysis  of,  this  method  see  Fix, 

Gunzburger  &  Nicolaides  [l3]»  Malkus  &  Olsen  [22],  and  Mercier  [24] 
and  references  therein. 

From  the  definition  II  it  follows  that  a  possible  method  to  achieve 
divergence-stability  of  relative  to  (if  it  is  not  already  there) 

is  to  increase  ,  or  decrease  (or  both).  The  well-known  Hood-Taylor 

[18]  element  is  obtained  from  continuous  piecewise  quadratic  velocities, 
piecewise  linear  pressures  exactly  this  way  by  requiring  the  pressures  also 
to  be  continuous.  The  element  analyzed  by  Arnold,  Brezzi  &  Fortin  [  1  ]  has  a  similar 
flavor:  for  velocities  the  continuous  piecewise  linears  are  enriched  by  cubic  bubble 

functions  (the  pressures  however  are  taken  to  be  continuous  piecewise  linear). 

In  the  context  of  the  equations  of  elasticity  it  has  been  proposed  to 
use  "reduced  integration"  on  the  term  involving  b  ,  while  maintaining 
simple  choices  for  both  the  meshes  and  Vh  .  We  refer  to  Zienkiewicz  [38],  Malkus 
&  Hughes  [21 ]  and  references  therein  for  computational  experience  with  reduced 
integration  methods,  and  an  account  of  their  relation  to  mixed  methods. 

An  obvious  question  at  this  point  seems  to  be: 

"are  the  problems  that  we  encountered  with  conforming 
piecewise  linear  elements  also  present  for  conforming 
piecewise  polynomials  of  degree  p  +  1  ,  p  j  1  ?" 

This  question  is  the  focus  of  the  next  two  sections. 
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6.  Divergence-stability  of  high-order  conforming  spaces 
Let  us  consider  the  spaces 

piecewige  polynomials  of  degree  $  p+1]^ 


(6.1) 


?„  -  [°p'p+11-V 

h  h 


{v  €  [ptp+ll.0]2  .  v  =  0  on  an } 

—  U  — 


on  an  arbitrary  family  of  triangulations  of  0  .  We  have  already  seen  that 

for  p  ■  0  these  spaces  are  not  in  general  divergence-stable.  Results  due  to 
de  Boor,  HtJllig  and  Jia  [11],  [19]  show  that  these  spaces  are  not  in  general 
divergence-stable  for  p  =  1  or  2  ;  this  lack  of  divergence-stability  is  well- 
documented  by  computational  experience  (cf.  section  7  for  more  details).  We 
now  turn  our  attention  to  the  case  p  £  3  .  In  order  to  state  rigorously  a 
result  concerning  the  divergence-stability  of  the  choice  (6.1)  for  p  >  3  we 
need  some  notation.  A  singular  internal  vertex  [25]  is  one  where  precisely 
four  triangles  meet  through  the  intersection  of  two  straight  lines,  as  shown 
in  Figure  3 .  A  singular  boundary  vertex  is  a  vertex  on  where  1  <  k  <  4 

triangles  meet  through  the  intersection  of  two  straight  lines.  There  are 


Figure  3 

Singular  internal  vertex, 
four  such  possibilities  as  shown  in  Figure  4. 


Figure  4 

Singular  boundary  vertices. 


Theorem  6.1  [ 32] 

Suppose  that  £  is  quasi-uniform  and  that  no  nonsingular  vertices 
degenerate  towards  singular  as  h  tends  to  zero.  Then  for  any  fixed  p  5  3 
the  spaces  (6.1)  are  divergence  stable  (in  the  sense  of  Definition  I). 
Furthermore  the  constant  c  entering  into  Definition  I  is  bounded  from  below  by 
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c'p  k  ,  c '  >  0 

with  c'  and  k  independent  of  h  and  p  . 

Remark  6.1 

By  quasi-uniform  we  mean  that  each  triangle  T  €  contains  a 
ball  of  radius  ph  ,  where  p  >  0  is  independent  of  h  .  Whether  this 
restriction  is  essential  to  guarantee  divergence-stability  of  the  spaces 
(6.1)  for  p  j  3  ,  we  do  not  know.  It  is  easy  to  see  that  the  spaces  (6.1) 
cannot  be  divergence-stable  if  a  nonsingular  vertex  tends  towards  becoming 

O 

singular  (as  seen  later  the  dimension  of  the  space  or  decreases 

by  one  when  a  vertex  becomes  singular).  This  does  not  imply  that  we  will 

observe  sub-optimal  convergence-rates  when  discretizing  the  Stokes 

problem  or  the  equations  of  elasticity,  v  -  1/2  ,  using  a  mesh  with  nearly-singular 

(nonsingular)  vertices,  however,  the  ratio  between  the  error  in  the  finite 

element  approximation  of  the  velocities  (or  displacements)  and  the  error 

of  the  best  (H^ -)  approximation  may  be  arbitrarily  large  depending  on  the  data.  As 

stated  in  Theorem  6.1  the  lower  bound  in  the  estimate  of  divergence-stability 

may  approach  zero  algebraically  in  p  ^  as  p  +  00  ,  we  do  not  know  whether 

this  is  indeed  the  case  or  whether  the  lower  bound  is  uniform  in  p  also. 

Note  that  divergence  stability  of  V,  requires  only  that  nonsingular  internal 

h  - 

vertices  do  not  degenerate  to  singularity  as  h  -*■  0  . 

For  the  proof  of  the  theorem,  as  well  as  more  details  concerning  the 
results,  we  refer  to  Scott  &  Vogelius  [32  ]  and  Vogelius  [37].  In  case  of 
the  equations  of  elasticity  (4.5)  Theorem  6.1  in  combination  with 
Proposition  5.4  leads  to  the  fact  that  the  spaces  (6.1)  give  quasi-optimal 
finite  element  approximations  to  (£V,  ^  ■  y~  V»_Uv)  for  fixed  p  *  3  .  The 
question  of  convergence  properties  thus  becomes  one  of  approximation  theory, 
and  in  addition  to  understanding  the  approximation  properties  of  it  is 
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important  to  understand  the  character  of  the  spaces  .  The  fact 

that  the  lower  bound  in  the  stability  estimate  approaches  zero  at  most 
algebraically  in  p  as  p  +  °°  was  used  in  Vogelius  [36]  to  prove  "almost" 
optimal  convergence  rates  for  the  so-called  p-version  of  the  finite  element 
method.  Using  an  interpolation  trick  it  was  shown  (on  a  smooth  domain,  with 
curved  elements  at  the  boundary)  that  the  p-version  converges  at  "almost"  optimal 
rates  in  the  energy  norm,  uniformly  with  respect  to  Poisson's  ratio.  "Almost" 
optimal  here  means  that  it  converges  at  any  rate  strictly  less  than  optimal.  A 
slight  variation  of  the  present  result  (permitting  for  curved  elements  at 
the  boundary)  could  likewise  be  used  to  verify  a  similar  result  for  combina¬ 
tions  of  the  h-  and  p-versions  (i.e.  when  simultaneously  changing  mesh  size  and 
degree  of  the  polynomials). 

For  the  case  of  the  Stokes  problem  it  is  clear  from  both  Proposition  5.2 
and  Proposition  5.3  that  in  addition  to  divergence-stability  and  good  approxi- 

O 

mation  properties  of  the  spaces  V,  it  is  important  to  have  good  approximation 

h 

properties  of  the  corresponding  pressure-spaces  11^  .  Since  Theorem  6.1 

O 

guarantees  divergence-stability  relative  to  any  subspace  of  V*V^  this 

O 

naturally  leads  to  the  question  of  finding  a  characterization  of  the  spaces  V*V,  . 

n 

Part  of  this  characterization  (its  necessity)  is  implicit  in  the  work  of  Nagtegaal, 
Parks  &  Rice  [27]  and  that  of  Mercier  [23]. 

Let  be  an  internal  vertex  where  four  triangles  meet  with  the  common 

edges  lying  on  either  the  x^-axis  or  the  x?-axis,  as  shown  in  Figure  5. 


Figure  5 


Let  v  €  and  set  =  vjy.  i  =■  1,...,4  .  Since  v  is  continuous 

14  i 

vanishes  identically  on  the  edge  e^  ,  consequently 


(6.2a) 


8  1  _9_  4 

9x^  V1  9x^  V1 


on  e^  . 


Similarly  we  get 


The  vertex  x^  is  common  to  all  the  edges  e^  through  e^  and  by  summation 
of  (6.2a-d)  at  Xg  we  8et 

(6.3)  7,Z1<*0)  +  7*^3(2So)  *  7*^2(^o)  +  7'Z^2o>  * 

which,  by  introduction  of  4>  ■  V*v  ,  may  be  restated 

4 

(6.4)  l  (-l)SlT  (O  -  0  . 

i=l  i 

The  chain  rule  of  differentiation  gives  that 

(6.5)  7*(T_1voT)(x)  -  (V*v) (Tx) 

2  2  -1 

for  any  invertible  affine  map  T  :  ]R  ->11  .  If  v  €  then  so  is  T  v(Tx) 

(on  the  corresponding  triangulation^ and  it  thus  follows  immediately  from  (6.5)  tha^ 
(6.3)  (or  (6.4))  must  be  satisfied  at  all  internal  vertices  which  can  be  obtained  by 
an  affine  transformation  of  one  as  in  Figure  5.  These  vertices  are  exactly  the 
singular  internal  vertices. 

This  simple  calculation  thus  shows  that  for  any  p  j  0 

(6.6)  7*7h  —  ^h^ *  ^  “  {(discontinuous)  piecewise 

polynomials,  $  ,  of  degree  s  p 
that  satisfy  (6.4)  at  singular 
internal  vertices}  . 

O 

For  the  spaces  ,  with  homogeneous  Dirichlet  boundary  conditions,  we  find 
two  additional  sets  of  contraints,  namely  with  <J>  *  V»v: 


(6.7) 


l  (-D  ♦  Ij-  (Xq)  =*  0 

i=l  i 


at  any  singular  boundary  vertex  x^  , 
where  1  <  k  <  4  triangles  meet  as  in 
Figure  4,  and 


(6.8) 


j  <j>  dx  =*  0 

J  o 


If  we  define 


p[p]>  1  xx  {(discontinuous)  piecewise  polynomials, 
h 

<p,  of  degree^  p,  that  satisfy  (6.4)  at  singular 
internal  vertices,  (6.7)  at  singular  boundary 
vertices  and  furthermore  satisfy  (6.8)} 


then  this  shows  that  for  any  p  >,  0 


(6.9) 


V*V  c  ptPl 
h  —  h 


It  was  shown  in  Vogelius  [37]  that  (6.6)  is  indeed  an  equality  provided  p  £  3  . 

By  a  similar  combinatorial  argument  it  is  proven  in  Scott  &  Vogelius  [32 ] 
that  (6.9)  is  also  an  equality  for  any  p  >,  3  (remember,  &  is  connected).  This 

O 

characterizes  the  spaces  V*V^  and  V»V^  on  an  arbitrary  triangulation  for  any  p  5  3. 
Remark  b.2 


From  the  characterization  nf  given  above  it  follows  that 


the  space 


(6.10) 


{continuous  piecewise  polynomials  of  degree  C  p 
that  satisfy  (6.8)} 
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is  a  subspace  of  V*V^  provided  p  >,  3  .  Selecting  the  pressure  space 
to  be  as  in  (6.10),  p  £  3,  generalizes  the  methods  studied  by  Hood  &  Taylor 
[18]  for  p  ■  1  and  2  .  Tneorem  6.1  shows  that  these  generatized  Hood- 
Taylor  elements,  P  5  3  ,  lead  to  quasi-optimal  approximation.  Bercovier 
&  Pironneau  [  3 ]  proved  the  same  to  be  the  case  for  the  Hood-Taylor  element 
for  p  *  1  .  It  would  be  natural  to  conjecture  that  the  case  p  *  2  also 
leads  to  quasi-optimal  approximations. 

Remark  6.3 

If  one  uses  the  approach  (5.5)  to  compute  the  finite  element  approxi¬ 
mation  to  the  velocity  IJ  in  the  Stokes  problem  then  it  is  important  to  have 
a  local  basis  for  Z^  .  Assuming  that  IT^  has  been  selected  so  that 
£  Z  it  then  follows  that,  for  all  p  J  0  , 

(6.11)  Z,  »  (v€V.  :  Vv  -  0}  *  V  x  p[P+2],l 

h  —  h  —  h 


where 


=  {C1  piecewise  polynomials  of  degree  < 

p  +  2  all  of  whose  first  derivatives 
vanish  at  the  boundary}. 


Note  that  elements  of  ptP+2],l  are  constant  on  each  connected  component 

h 

of  3ft  ,  however  the  constants  need  not  be  the  same;  if  3ft  is  connected  then  we 
shall  always  pick  the  constant  to  be  zero  for  functions  in  , 

A  basis  for  P^p+^,\  p  i  3  ,  (with  no  boundary  conditions)  was  constructed 
in  Morgan  &  Scott  [25];  the  corresponding  basis  for  p^P+^»^  (  p  j  3  ,  is  not 
altogether  obvious,  but  is  described  in  Scott  &  Vogelius  [32],  The  finite 

O  O 

element  method  corresponding  to  the  spaces  V,  and  IT,  =  V*V,  in  the  case 

n  n  n 

p  ■  3  is  currently  being  tested  numerically  for  the  full  (nonlinear,  time- 
dependent)  Navier-Stokes  equations. 


Results  concerning  lower-degree  spaces 


There  are  really  three  aspects  of  the  finite  element  approximation  to 
the  Stokes-equations  or  the  equations  of  elasticity  addressed  in  the  previous 
sections . 

a)  Bounds  for  a  maximal  right  inverse  for  the  divergence 

O 

operator  on  and  ,  or  estimates  for  the 

divergence-stability  of  these  spaces  relative  to 
families  of  pressure  spaces  11^  . 

O 

b)  Characterization  of  the  spaces  V*V^  or  . 

c)  Determination  of  the  approximation  properties  of  the 
spaces  Vh  fl  {v:V*v=0>  (or  fi  {v:V*v=0}). 

For  the  equations  of  elasticity  (4.5)  we  are  as  indicated  by 
Proposition  5.4  mainly  concerned  with  a)  and  b)  whereas  the  aspect  c)  is  also 
of  importance  when  solving  the  Stokes-problem  (4.4).  As  pointed  out  in 
sectioiB  5  and  6  (particularly  in  Remark  5.2  and  Remark  6.1)  c)  may  be  of  independent 
interest  even  though  a)  and  b)  do  not  have  satisfactory  answers.  We  shall  now 
discuss  the  aspects  a)-c)  for  the  spaces 


(7.1) 


V.  =  [ptP+1J’°]2  and  V,  -  [plP*1).0]2  § 


P  *  0,1,2  . 


We  know  of  very  few  results  that  are  valid  on  a  quite  arbitrary  triangulation 

for  these  low  degree  spaces;  as  a  consequence  we  shall  restrict  our  attention 

to  the  triangulations,  a  local  picture  of  which  are  shown  in  Figure  1  and 

X  2 

Figure  2  (we  denote  these  and  £  respectively).  On  these  triangula- 

tlons  the  dimension  formula  conjectured  by  Strang  [33] 

(7.2)  dim(P^P+2]’1)  -  y(p+3)(p+4)T-(2p+5)E0+3V0+a0 
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is  known  to  hold  also  for  p  *  0,1  and  2  (cf.  Morgan  &  Scott  [26]).  Here 

T  denotes  the  number  of  triangles,  Eg  is  number  of  internal  edges,  Vg 

is  the  number  of  internal  vertices  and  Og  is  the  number  of  singular 

internal  vertices.  On  a  simply  connected  domain  the  null  space  of  the 

divergence  operator  acting  on  V,  ■  t^e  curi  Qf  ptp+2]»l  . 

n  h  h 

it  thus  follows  as  in  Vogelius  [37]  (or  Scott  &  Vogelius  [32])  that  whenever 

(7.2)  holds  7»V^  must  have  the  same  dimension  y(p+2) (p+l)T  -  Og  as  the 
space  p|p^’  1  (cf.  (6.6)).  Since  7*V,  c  P^P^’  ^  we  conclude  that 

(7.3)  7*V,  =  P/Pl*-1 

h  h 

whenever  (7.2)  holds  for  a  given  triangulation,  and  in  particular  for  p  *  0,1 

rl  r2 

and  2  on  the  triangulations  and  4^  •  (For  P  “  2  we  could  indeed 

have  derived  (7.2)  for  much  more  general  triangulations,  cf.  Morgan  &  Scott 

[26].)  By  a  hole-filling  procedure  we  can  extend  our  argument  to  verify 

(7.3)  even  though  the  domain  is  not  simply  connected.  The  characterization 

2 

(7.3)  for  the  case  p  ■  0  and  the  triangulation  ^  was  also  noted  by  Fix, 
Gunzburger  &  Nicolaides  (cf.  [13]). 

In  contrast  the  relation  between  the  spaces  7*V^  and  P^1^’  ^  is  not 
nearly  as  simple  for  p  •  0,1  and  2  as  the  characterization  given  in  the 
previous  section,  for  p  %  3  ,  might  lead  one  to  believe.  Let  us  start  by 
considering  the  piecewise  linear  case,  p  ■  0  .  If  o  denotes  the  total  number 

of  singular  vertices,  including  singular  boundary  vertices,  then  we  have  for 

rl  r2 

the  triangulations  ^  and  ^ 

dim  Vh  -  2Vq  and  dim  -  T-o-1 

(excluding  the  trivial  case  of  a  rectangle  divided  into  two  triangles  by  the 
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0  -roi  -i 

diagonal,  when  both  V,  and  P£  *  consist  of  0  only).  Using  the  relations 


E-Eq  -  V-VQ  ,  E+Eq  =  3T 


and  Euler's  formula 


T+V-E  =  1 


(assuming  ft  is  simply  connected),  we  get 


(7.4) 


dim  p/03’”1  -  dim  V 
h  h 


(E-Eq)  -  o  -  3  ; 


T  ,  Eq  ,  Vq  are  as  before  and  E  denotes  the  total  number  of  edges,  V  the 

total  number  of  vertices.  From  (7.4)  we  immediately  conclude  that  if 

0  ,~f01  -1  ° 

o  <  E  -  Eq  -  3  then  dim  <  dim  and  consequently  is  a 

proper  subspace  of  ^  (see  Nagtegaal,  Parks  &  Rice  [27]  and  Malkus  & 

Hughes  [21]  for  a  similar  constraint-counting  method).  If  ft  is  a  rectangle 
cl 

and  ^  is  used  for  a  triangulation  then  there  are  exactly  2  singular  vertices 

0 

(one  in  the  upper-left  and  one  in  the  lower-right  corner  of  ft  );  thus  V*V, 

h 

is  a  proper  subspace  of  ^ ^  except  for  the  trivial  case  that  ft  is 

O  O 

divided  into  only  two  triangles.  Since  V,  fl  Z  *  {v€V,  :  V*v=0}  **  {0}  on  the 

n  —  n  — 

rl  0 

mesh  it  follows  that  7*  is  injective  on  and  hence  we  get  from 


(7.4)  that 


dim  Pi01’-1  -  dim  7*V.  -  (E-En)-5 
h  h  0 


on  a  rectangle  triangulated  by  the  mesh  ^  (excluding  the  trivial  case  of  only 


2  triangles) . 


In  the  case  that  ft  is  a  rectangle  divided  into  triangles  by  the 

mesh  7?  then  the  dimension  of  p/^*^  *  p/^’^  fl  ■  0  on  3ft} 

‘‘ti  h  h  on 


.SVV-  , 


‘  •/  %’  *.  »/-  <i  •  .  •  »  •  . 


.  **•  •*»  **,  »*,  <*, »' 


has  been  calculated  by  Chui,  Schumaker  &  Wang  [ 8  ]  to  be 


dim  p£p1,1  -  <x  -  (E-EQ)  +  4  , 

provided  Q  has  at  least  two  boundary  edges  on  each  side.  As  a  consequence 
of  this  and  (6.11)  it  follows  that 

dim  Vu  -  dim  7-V.  +  a  -  (E-E  )  +  4  , 
n  n  u 

and  thus,  using  (7.4),  we  find  that 

(7.5)  dim  PJ01*'1  -  dim  7»V.  -  1 

n  n 

r2 

for  the  choice  (7.1),  with  p  -  0  ,  on  this  mesh  ^  . 

Now  consider  the  case  of  piecewise  quadratic  fields  ,  i.e.  the  choice 
Vh  ■  (pj^’^)2  on  the  mesh  ^  .  The  dimension  of  P^^ :  iji  = 

3H) 

—  "0  on  30}  on  a  rectangle  triangulated  by  this  mesh  has  been  determined 
in  Chui,  Schumaker  &  Jang  [  9  j,  and  via  calculations  simular  to  those  just 
given  we  find  that 

(7.6)  dim  P/11’"1  -  dim  7*V.  -  3  . 

n  n 

Note  that  this  agrees  with  the  results  of  Malkus  &  Olsen  [22  ] .  For  the  space 

•  oral  A  a  a 

Vh  -  (P^  '*  )  on  £  the  situation  is  different.  Mercier  [23]  observed 
that  in  this  case  has  a  local  basis  given  by  the  Fraeij  s  de  Veubeke- 

Sander  cubic  macro-element.  Using  calculations  similar  to  those  above  one 
can  show  that 

(7.7)  7-V.  -  P*11’-1 

n  n 

2 

on  a  rectangle  triangulated  by  the  mesh  (p«l). 

For  the  case  p  »  2  ,  that  is  piecewise 
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0  0  f  3  I  0  2 

cubic  fields,  Vh  *  (Pl  *  )  ,  with  homogeneous  boundary  conditions  on  a 

rectangle  triangulated  by  the  mesh  L  %one  can  use  the  local  basis  for  ?l4J* 

n  *  h 

given  via  macro-elements  in  Douglas  et  al.  [12],  to  show  that 


(7.8) 


7*7, 


p[2],-l 

h 


To  the  best  of  our  knowledge,  no  result  is  known  for  the  mesh  for  p  =  2  , 

although  it  seems  clear  that  one  would  not  expect  (7.8)  to  hold  in  general. 

Whenever,  as  in  (7.5)  and  (7.6),  the  codimension  of  7»V,  in  ^ 

n  h 

is  non-zero  but  independent  of  the  mesh  size,  it  is  reasonable  to  conjecture 
that  the  orthogonal  complements  consist  of  global  modes.  For  an  explicit 

O 

calculation  of  the  global  constraint  on  V in  the  case  corresponding  to 
(7.5),  see  Olsen  [28]. 

The  previous  discussion  centered  on  the  characterization  of  7*V,  or 

n 

O 

7*V^  for  low  degree  polynomials;  we  now  turn  to  the  aspects  a)  and  c)  listed 
at  the  beginning  of  this  section.  Proposition  5.1  shows  that  there  is  a  close 

O 

connection  between  approximation  properties  of  the  spaces  V  D  {v  :  7*v=0}  (or 
Vh  n  W  :  7*v=0} )  and  the  divergence  stability  of  Vh  (or  Vh  ,  respectively). 

In  the  following  we  shall  for  simplicity  assume  that  ft  is  simply 

connected.  Then  V,  D  {v  :  7*v»0)  ■  7  x  ptp+2l»l  and  v  fl  (v  :  7*v  *  0} 

h  —  —  h  h  —  — 

*  7  x  f  and  this  in  combination  with  Proposition  5.1  shows  that  there  is  a 

close  connection  between  approximation  properties  of  an(j  ptp+2l»l  and 

h  h 

divergence-stability.  The  next  proposition  elaborates  more  on  that  connection. 
Proposition  7.1 

Let  p  *  0,1  or  2  and  assume  that  there  exist  €  Hp+3(ft)  , 
c  >  0  ,  and  0  <  ot  <  p  +  3  such  that 


L2(ft)^ 


c  hP+3'a 


inf  1 1 


where  Che  inf  is  taken  over  <*>  €  P^H  ''  (for  some  quasi-uniform  family  of 

h 

triangulations).  Then 

b(v,q)  g 

inf  sup  -  <  Ch 

q€V-V  \{0}  v*V  1(0}  ||v||  ,  ||q|lL 
h  h  Hl(£j) 


for  the  subspaces  «  (P^^^*^)2  ,  with  8 
Proof : 

By  contradiction  let  us  assume  that 


p+1 

P+3  ° 


h-6  inf 


b(v,q) 


,ev.vhvto)  v«h\{o> 


is  unbounded  as  h  -*•  0  .  Using  the  same  argument  as  in  the  proof  of  Proposition  5.1 


we  get  directly  from  the  definition  of  that 


(7.9) 


^€V^ll{v:V»v*0} 


| j U— z  1 1  1  <  CC'V6  inf  ||u-v|| 


h  (ft) 


v€V. 

—  n 


H  (ft) 


for  any  U  €  [H1  (ft) ] 2  0  {v  :  V*v  -  0}  .  For  any  <l>  €  Hp+3(ft)  we  get  that 
7  X  |(|  (  [HP+2(J1)]2  n  (v  :  V *v  =>  0},  and  by  insertion  into  (7.9)  it  follows 


inf  ||vxip-Vx$||  $  CC~1h"8  inf  ||Vxi|>-v||  . 

^p[p+2],l  HX(ft)  h  v€Vh  H  (P.) 

.  n  VP+1  r_1>rB  11,1.  II 


*  C  h^  chxh  p  jk!l  . 

HP  J(ft) 


From  this  we  conclude  that,  for  any  €  HP  3 (ft)  , 


,-l,  p+!-B  i|,  |i 


’ll 


Using  the  results  of  Bramble  and  Scott  [5]  (these  results  are  valid  if  we  require 
that  all  corners  of  ft  have  interior  angles  <  2  it  )  we  thus  get 


,tJ«l.il,^llL,<“>  *  IM!^ 


which  due  to  the  fact  that  is  unbounded  as  h  0  produces  a  contradiction 

to  the  assumptions  of  this  proposition. 

Proposition  7.1  could  just  as  easily  have  been  verified  for 

HP+3(ft)flH2(ft)  ,  p';P+2]’1  =  {$epfP+2]*1  ;  *  =  |±  =  o  at  3ft} 
n  n  3n 

and  V.  -  (ptpd-H.0,2 
n  n 

In  the  light  of  Propositions  5.1  and  7.1  we  shall  center  our  discussion 
of  the  aspects  a)  and  c)  on  known  results  concerning  the  approximation  rates 
of  P^p+2]’1  and  P^p+2J'1  on  the  meshes  and  £2  . 

First,  suppose  p  *  0  ,  the  case  of  piecewise  linear  fields.  On  the 
mesh  we  know  as  previously  stated  that 

?l2U1  =  (0}  ; 

n 

this  space  obviously  has  no  approximation  rate  associated  to  it  and  the 
assumptions  of  Proposition  7.1  are  thus  satisfied  with  a  =  3  ,  i.e.  the 
"best"  constant  in  the  "divergence-stability"  estimate  is  s  C  h  ,  cf.  Gunzburger 
&  Nicolaides  [16].  We  do  not  anticipate  that  relaxing  the  boundary  conditions 
would  lead  to  any  approximation  rates  for  the  spaces  P^2^*3  on  the  mesh  ^  , 
though  we  do  not  know  this  for  a  fact.  On  the  mesh  ^  Fix,  Gunzburger  & 
Nicolaides  [13]  give  a  construction  of  a  uniformly  bounded  inverse  for  the 
divergence  operator  P^  *”1  +  \  m  [P^’0]2  »  uniformly  bounded  that 

is  in  8(H~^;L_)  ,  not  8(L2;H1')  as  we  are  concerned  with.  Powell  [ 30 ]  and 


.V.V.V.V /.-V". 


mown 


Mercier  [23]  both  show  that  P^  ’  has  optimal  approximation  rates  on  this 

mesh.  Malkus  and  Olsen  [22]  coni ecture  that  optimal  approximation  properties 
o  r  21  1  2 

hold  as  well  for  Pl  ’  on  ^  ,  although  this  has  not  been  proven. 

”[211 

(Because  of  the  local  basis  for  P  *  given  by  Powell  [30],  one  would  have 
a  *  1  at  worst.) 

Now  consider  piecewise  quadratic  fields,  i. e.  p  =  1  .  On  the  mesh 
de  Boor  and  Hdllig  [11]  show  that  approximation  by  is  suboptimal 

by  precisely  one  order  of  h  in  This  result  (adapted  to  L^)  accord¬ 

ing  to  Proposition  7.1  leads  to  the  conclusion  that  the  "best"  constant  in  the 

"divergence-stability"  estimate  for  the  piecewise  quadratic  spaces  is 

1/2 

$  C  h  .  Malkus  and  Olsen  [22]  report  numerical  evidence  that  the  corresponding 

O 

constant  for  the  piecewise  quadratic  spaces  (with  boundary  conditions)  is 

0(h)  on  the  mesh  .  On  the  mesh  ^  the  situation  is  quite 

[3 ]  2  *r31  i 

different  since  both  P^  ’  and  P^  J ’  are  long  known  to  have  optimal 

approximation  rates,  cf.  Mercier  [23];  indeed,  the  recent  results  of 

Boland  and  Nicolaides  [4  ]  together  with  the  result  (7.7)  show  that  the 

« 

corresponding  spaces  and  are  divergence  stable  on  the  mesh  for 

p  =*  1  .  (For  macroelement,  one  takes  simply  the  quadrilateral  surrounding 
each  singular  vertex.  Local  stability  is  then  guaranteed  by  (7.7)  applied  to 
each  macroelement.  The  reason  that  these  macroelements  are  locally  stable 
in  the  sense  of  Boland  and  Nicolaides  [  4 ]  is  that  the  boundary  vertices  of 
each  macroelement  are  nonsingular,) 

_1 

For  piecewise  cubic  fields  (p  ■  2)  on  the  mesh  ,  numerical  exper- 

O 

ience  (e.g.  Szabo  et_  a]^.  [34])  had  indicated  that  neither  nor  was 

divergence-stable.  Moreover,  the  results  of  Jia  [19],  similar  to  those  of 

de  Boor  and  Hollig  [11]  quoted  above,  also  lead  to  this  conclusion  via 

-2 

Proposition  7.1.  For  the  mesh  and  p  =  2  ,  the  results  of  Douglas,  et_  al .  [12  ] 
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prove  that  (7.8)  holds  and  consequently  that  both  and 

O 

V,  are  divergence-stable,  in  the  same  fashion  as  for  the  case  p  =  1 

described  previously.  Finally,  note  that  all  the  results  for  p  =  1  and 
c2 

2  for  the  mesh  ^  hold  as  well  for  a  mesh  based  on  the  macroelement  of 
Clough  and  Tocher,  c£.  Mercier  [23], 

To  summarize  the  previous  discussion: 

1)  The  characterization  of  the  range  of  the  divergence  operator 

on  spaces  of  continuous  piecewise  polynomials  of  degree  at  most  p  +  1  ,  p  >,  3  , 
given  in  section  6,  is  widely  valid  also  for  p  £  2  ,  provided  no  boundary 
conditions  are  imposed.  With  boundary  conditions  imposed  this  characterization 
fails  on  the  most  natural  triangulations. 

2)  The  divergence  stability  of  the  spaces  [p^P+^ ^  or  [f^P 

is  intimately  connected  to  the  approximation  properties 

of  the  spaces  pfp+2]»l  or  ptp+2],l  __  Since  these  spaces  of  piecewise 

h  h 

polynomials  have  essentially  one  additional  degree  of  freedom  for  each  singular 
vertex,  it  is  only  natural  that  the  spaces  [p^P+-*-^  or  [ptP+l]’0j2  }  p  <:  2  ,  are 

much  more  likely  to  be  divergence-stable,  the  more  singular  vertices  the 
triangulation  has.  In  this  sense  singular  vertices  are  desirable  when  working 
with  piecewise  polynomials  of  degree  s  3. 


References 


[1]  Arnold,  D.  N. ,  Brezzi,  F.  and  Fortin,  M.  ,  "A  stable  finite  element  for  the 
Stokes  equations".  Calcolo,  to  appear. 

[2]  BabuSka,  I.,  "Error-bounds  for  the  finite  element  method",  Numer .  Math, 
vol.  16  (1971),  pp.  322-333. 

[3]  Bercovier,  M.  and  Pironneau,  0.,  "Error  estimates  for  the  finite  element 
method  solution  of  the  Stokes  problem  in  the  primitive  variables", 

Numer .  Math,  vol.  ^33  (1979),  pp.  211-224. 

[4]  Boland,  J.  M.  and  Nicolaides,  R.  A.,  "Stability  of  finite  elements  under 
divergence  constraints",  SIAM  J.  Numer.  Anal .  vol.  20  (1983),  pp.  722-731. 

[5]  Bramble,  J.  H.  and  Scott,  R. ,  "Simultaneous  approximation  in  scales  of 
Banach  spaces",  Math.  Comp .  vol.  32_  (1978),  pp.  947-954. 

[6]  Brezzi,  F. ,  "On  the  existence,  uniqueness  and  approximation  of  saddle- 
point  problems  arising  from  Lagrangian  multipliers",  R.A.I.R.O.  vol.  8^ 

R2  (1974),  pp.  129-151. 

[7]  Chui,  C.  K.  and  Wang,  R.-H. ,  "On  smooth  multivariate  spline  functions". 
Math.  Comp . ,  to  appear. 

[8]  Chui,  C.  K. ,  Schumaker,  L.  L.  and  Wang,  R.  H. ,  "On  spaces  of  piecewise 
polynomials  with  boundary  conditions,  II.  Type-1  triangulations",  to 
appear  in  Proceedings  of  the  second  Edmonton  Conference  on  Approximation 
Theory. 

[9]  Chui,  C.  K. ,  Schumaker,  L.  L.  and  Wang,  R.  H. ,  "On  spaces  of  piecewise 
polynomials  with  boundary  conditions.  III.  Type-2  triangulations",  to 
appear  in  Proceedings  of  the  second  Edmonton  Conference  on  Approximation 
Theory. 

[10]  Crouzeix,  M.  and  Raviart,  P.-A.  ,  "Conforming  and  nonconforming 
finite  element  methods  for  solving  the  stationary  Stokes  equations", 
R.A.I.R.O.  vol.  1_  R3  (1973),  pp.  33-76. 


[11]  de  Boor,  C.  and  HiJllig,  K. ,  "Approximation  order  from  bivariate  C^-cubics 
A  counterexample",  MRC  Tech.  Rep.  // 2389  (1982),  to  appear  in  Proc . 

Amer .  Math.  Soc . 

[12]  Douglas,  Jr.,  J.,  Dupont,  T.,  Percell,  P.,  and  Scott,  R. ,  "A  family  of 

finite  elements  with  optimal  approximation  properties  for  various 
Galerkin  methods  for  2nd  and  4th  order  problems",  R.A.I.R.O.  Num.  Anal . 
vol •  13  (1979),  pp.  227-255. 

[13]  Fix,  G.  J. ,  Gunzburger,  M.  D.,  and  Nicolaides,  R.  A.,  "On  mixed  finite 
element  methods  for  first-order  elliptic  systems",  Numer.  Math,  vol. 
(1981),  pp.  29-48. 

[14]  Glowinski,  R. ,  Numerical  Methods  for  Nonlinear  Variational  Problems, 

2nd  Ed.,  Springer-Verlag  Series  in  Comp.  Physics,  to  appear. 

[15]  Grisvard,  P. ,  Boundary  value  problems  in  non-smooth  domains.  Part  2. 
University  de  Nice,  France,  1981. 

[16]  Gunzburger,  M.  D.  and  Nicolaides,  R.  A.,  "The  computational  accuracy  of 
some  finite  element  methods  for  incompressible  viscous  flow  problems", 
ICASE  Report  No.  81-7. 

[17]  Gunzburger,  M.  D.  and  Peterson,  J.  S.,  "On  conforming  finite  element 
methods  for  the  inhomogeneous  stationary  Navier-Stokes  equations", 

Numer.  Math,  vol.  42  (1983),  pp.  173-194. 

[18]  Hood,  P.  and  Taylor,  C.,  "Navier-Stokes  equations  via  mixed  interpolation 
in  Finite  Elements  in  Flow  Problems,  Oden  et  al. ,  eds. , 

Univ.  of  Alabama  Press,  1974,  pp.  121-132. 

[19]  Jia,  R. ,  "Approximation  by  smooth  bivariate  splines  on  a  three-direction 
mesh",  MRC  Tech.  Rep.  #2494  (1983). 

[20]  Kellogg,  R.  B.  and  Osborn,  J.  E.,  "A  regularity  result  for  the  Stokes 
problem  in  a  convex  polygon".  Journal  of  Functional  Analysis ,  vol.  21_ 

(1976),  pp.  397-431. 


[21]  Malkus,  D.  S.  and  Hughes,  T.  J.  R. ,  "Mixed  finite  element  methods- 
reduced  and  selective  integration  techniques:  a  unification  of  concepts 
Comp.  Meth.  Appl .  Mech.  Eng. ,  vol.  15.  (1978),  pp.  63-81. 

[22]  Malkus,  D.  S.  and  Olsen,  E.  T. ,  "Obtaining  error  estimates  for  optimally 
constrained  incompressible  finite  elements".  Comp.  Meth.  Appl .  Mech. 

Eng. ,  to  appear. 

[23]  Mercier,  B.,  "A  conforming  finite  element  method  for  two  dimensional, 
incompressible  elasticity",  Int .  J_.  Numer .  Meths .  Eng,  vol.  14^  (1979), 
pp.  942-945. 

[24]  Mercier,  B.,  Lectures  on  Topics  in  Finite  Element  Solution  of  Elliptic 
Problems,  Springer-Verlag,  1979. 

[25]  Morgan,  J.  and  Scott,  R. ,  "A  nodal  basis  for  C1  piecewise  polynomials 
of  degree  n  *  5",  Math.  Comp,  vol.  29^,  (1975),  pp.  736-740. 

[26]  Morgan,  J.  and  Scott,  R. ,  "The  dimension  of  the  space  of  C^  piecewise 
polynomials",  unpublished  manuscript,  1975. 

[27]  Nagtegaal,  J.  C.,  Parks,  D.  M.  and  Rice,  J.  R. ,  "On  numerically  accurate 
finite  element  solutions  in  the  fully  plastic  range",  Comp.  Meth.  Appl . 
Mech.  Eng.,  vol.  4.  (1974),  pp.  153-177. 

[28]  Olsen,  E.  T. ,  "Stable  finite  elements  for  non-Newtonian  flows:  first- 
order  elements  which  fail  the  LBB  condition",  Ph.D.  Thesis,  Illinois 
Inst,  of  Tech.,  1983. 

[29]  Osborn,  J.  E.,  "Regularity  of  solutions  of  the  Stokes  problem  in  a  poly¬ 
gonal  domain".  In  Numerical  Solution  of  Partial  Differential  Equations, 
III,  B.  Hubbard  ed..  Academic  Press,  1976,  pp.  39j-411. 

[30]  Powell,  M.  J.  D. ,  "Piecewise  quadratic  surface  fitting  for  contour 
plotting",  in  Software  for  Numerical  Mathematics,  D.  J.  Evans,  ed.. 
Academic  Press,  1974,  pp.  253-271. 


[31]  Pritchard,  W.  G. ,  Renardy,  Y.,  and  Scott,  L.  R. ,  "Tests  of  a  numerical 
method  for  viscous,  incompressible  flow.  I:  Fixed-domain  problems", 
to  appear. 

[32]  Scott,  L.  R.  and  Vogelius,  M. ,  "Norm  estimates  for  a  maximal  right 
inverse  of  the  divergence  operator  in  spaces  of  piecewise  polynomials". 
Tech.  Note  BN-1013,  Inst.  Phys.  Sci.  &  Tech.,  Univ.  of  Maryland,  1983. 

[33]  Strang,  G. ,  "Piecewise  polynomials  and  the  finite  element  method". 

Bull.  AMS  vol.  79  (1973),  pp.  1128-1137. 

[34]  Szabo,  B.  A.,  Basu,  P.  K. ,  Dunavant,  D.  A.  and  Vasilopoulos,  D., 

Adaptive  finite  element  technology  in  integrated  design  and  analysis. 
Report  WU/CCM-81/1,  Washington  University,  St.  Louis. 

[35]  Temam,  R. ,  Navier-Stokes  Equations,  North-Holland ,  1977. 

[36]  Vogelius,  M. ,  "An  analysis  of  the  p-version  of  the  finite  element  method 
for  nearly  incompressible  materials.  Uniformly  valid,  optimal  error 
estimates",  Numer.  Math,  vol.  4_1_  (1983),  pp.  39-53. 

[37]  Vogelius,  M. ,  "A  right-inverse  for  the  divergence  operator  in  spaces  of 
piecewise  polynomials.  Application  to  the  p-version  of  the  finite  element 
method",  Numer.  Math,  vol.  41^  (1983),  pp.  19-37. 

[38]  Zienkiewicz,  0.  C.,  The  Finite  Element  Method  in  Engineering  Science, 
McGraw-Hill,  1971. 


